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In a quasi two-dimensional electron system with non-zero layer thickness, a parallel magnetic 
field can couple to the out-of-plane electron motion and lead to a severe distortion and eventual 
splitting of the Fermi contour. Here we directly and quantitatively probe this evolution through 
commensurability and Shubnikov-de Haas measurements on electrons confined to a 40-nm-wide 
GaAs (001) quantum well. We are able to observe the Fermi contour splitting phenomenon, in good 
agreement with the results of semi-classical calculations. Experimentally we also observe intriguing 
features, suggesting magnetic-breakdown-type behavior when the Fermi contour splits. 


In a strictly two-dimensional electron system (2DES) 
with zero layer thickness, the electron in-plane motion is 
unaffected by a parallel magnetic field (S| |) ■ However, for 
a quasi-2DES, such as electrons in a quantum well (QW) 
with finite width, can couple to electrons’ out-of-plane 
motion, thus also affecting their in-plane motion. This 
can have profound consequences. For example, the Fermi 
contour, which is circular in an isotropic system such as 
the 2DES in GaAs QWs, becomes severely distorted by 
and could even split into two, tear-drop shaped con¬ 
tours if B\\ is sufficiently strong (Fig. 1). B\\ also causes 
spin-polarization, leading to the formation of two distinct 
Fermi contours with different enclosed areas for different 
spins. Therefore, the spin-degenerate Fermi contour at 
B = 0 could split into two pairs of smaller contours in 
the presence of a large i?|| (Fig. 1). Since the shape of 
Fermi contour, when rotated by 90°, reflects that of the 
cyclotron orbit in real space [T], its evolution as a func¬ 
tion of iJ|| corresponds to the evolution of the electron 
trajectory. Understanding this i3||-induced Fermi con¬ 
tour splitting is of fundamental importance, specially for 
spintronic devices where application of iJ|| is often used 
for spin-polarization nisi- 

Several transport studies on 2DESs, confined to cou¬ 
pled double- [m and triple-QW systems [10] have pre¬ 
viously explored the Fermi contour splitting. In these 
studies, features such as kinks in the H||-induced mag¬ 
netoresistance and inter-layer tunneling were associated 
with the splitting. Electrons in very wide single QWs, 
which are essentially bilayer systems, also produced sim¬ 
ilar results m- Compared to the earlier works, our study 
here incorporates the following novelties: (1) The 2DES 
is confined to a single QW with a single-layer-like charge 
distribution at B\\ = 0 (see inset of Fig. 1(A)). (2) We 
probe the splitting of the Fermi contour for both spin 
species via Shubnikov-de Haas (SdH) oscillations. (3) 
We use measurements of commensurability oscillations 
(COs), also known as Weiss oscillations [H], to directly 
map out the Fermi contour and capture its distortion and 
the eventual splitting. 

Our sample, grown via molecular beam epitaxy, is a 40- 


nm-wide, GaAs (001) QW which is located 190 nm under 
the surface. The QW is flanked on each side by 95-nm- 
thick Alo.24Gao.76As spacer layers and Si 5-doped layers. 
The 2DES density is n = 1.75 x 10^^ cm“^, and the 
mobility is ~ 20x10® cm^/Vs. We fabricated a strain- 
inducing superlattice with a period a = 400 nm on the 
surface of our sample, an L-shaped Hall bar (Fig. 2 (a)). 
The superlattice, made of negative electron-beam resist, 
modulates the potential through the piezoelectric effect 
in GaAs [T5HT7] . For iJ||-dependent measurements, we 
first apply a large B-field in the plane along [110]. The 
sample, mounted on a single-axis tilting stage, is then 
slowly rotated around [110] using a computer-controlled, 
brushless DC motor to introduce a small component of 
the field perpendicular to the 2D plane. This B± induces 
SdH oscillations in the unpatterned reference and COs 
in the modulated regions of the Hall bar. (Note that 
iJ|l ~ B because ^ B.) We pass current along the L- 
shaped Hall bar and measure the longitudinal resistances 
simultaneously for both arms. B± is extracted from a 
linear fit of the Hall resistance measured in the reference 
region. All measurements are carried out at 300 mK. 

Figures 1(A)-(E) highlight the key points of our study. 
The Fermi contours are derived from calculations based 
on an 8 X 8 Kane Hamiltonian with no adjustable pa¬ 
rameters m- We include B|| = {Bx,By,Q) via the 
vector potential A{z) — {zBy,—zBx,0) so that the in¬ 
plane canonical momentum k = (kx,ky,0) remains a 
good quantum number. The occupied states enclosed 
by the Fermi contours define the charge density from 
which the Hartree potential is derived self-consistently. 
As seen in Fig. 1(A), the Fermi contours of the two spin- 
subbands, which are identical for H||=0, start out circular 
m- As iJ|| increases, the 2DES becomes progressively 
more spin-polarized, thus two distinct spin Fermi con¬ 
tours are formed. The application of i?|| along [110], also 
elongates and splits the contours in the [110] direction. 
Figures 1(G) and 1(D) illustrate the splitting of the mi¬ 
nority and majority spin contours, respectively. The split 
contours separate even further along [110] under stronger 
H|| (Fig. 1(E)). The insets of Figs. 1(A)-(E) show how 
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FIG. 1. (color online) (A)-(E) Calculated Fermi contours for a 2DES with density n = 1.75 x 10^^ cm ^ confined to a 
40-nm-wide GaAs QW. B|| is applied along [110]. Insets show the corresponding charge distributions in light blue. Inset of (E) 
also shows the charge distribution from fc[110] > 0 (shaded pink region) and fc[110] < 0 (shaded yellow region) states. 




FIG. 2. (color online) (a) Sample schematics. The electron- 
beam resist grating covering the surface of each Hall bar arm 
is shown as blue stripes. Part of the [110] arm is left unpat¬ 
terned as a reference region, (b) B||-induced magnetoresis¬ 
tance for Rref- (A)-(E) mark the values that correspond 
approximately to the calculations of Fig. 1. 


the charge distribution along the confinement direction 
gradually evolves from a single-layer into a bilayer, with 
each of the “layers” corresponding to one part of the split 
Fermi contour (see Fig. 1(E) inset). 

As an indication of Fermi contour splitting, we first 
present the B||-magnetoresistance trace from the refer¬ 
ence region of the Hall bar (Fig. 2(b)). The trace demon¬ 
strates non-monotonic transport behavior. In particular, 
there are two pronounced kinks at = 11 T and = 14 T, 
marked by red and green circles, respectively. The po¬ 
sitions of these kinks agree well with the expected split¬ 
ting of the minority and majority spin contours (see Figs. 
1(C) and (D)). A kink in the i?||-magnetoresistance has 
been previously associated with the splitting of Fermi 
contour [5]. Here, we observe two kinks, suggesting a 
spin-dependent splitting of Fermi contours. 

We further investigate the splitting via SdH oscillations 
which directly probe the area enclosed by the Fermi con¬ 
tour. We expect that the splitting would be reflected 
as a jump in the SdH frequency. Figure 3(a) shows the 
SdH oscillations at different H|| while the corresponding 
Fourier transforms (FTs) are shown in Fig. 3(b). For 
H|| = 0 T, we observe two peaks, the stronger of which 
if SdH — 3-56 T) is for the spin-unresolved SdH oscilla¬ 
tions (marked by a dotted vertical gray line). The weaker 


peak at 7.08 T is very close to the value of 2/^^^ (marked 
by a solid vertical gray line) and corresponds to the spin- 
resolved oscillations [20]. Around H|| = 9 T, the spin- 
unresolved peak splits, with the lower frequency peak 
f~ (red square) corresponding to the electron density of 
the minority-spin-subband and the higher frequency peak 
/■*■ (green square) to the majority-spin-subband. Then, 
starting at i?|| = 11 T, another low-frequency peak 
(red square) appears at approximately /“/2, signaling 
the splitting of the minority-spin contour. The peak 
remains dominant between i3|| = 11.5 and 13.5 T where 
both /■'■ and f~ become very weak and essentially van¬ 
ish. However, at H|| = 15 T, another peak (green 
square) appears to the right of /^2 becomes the 
dominant feature in the FT spectrum up to H|| = 18 T. 
The sum of dose to implying that 

ft /2 originates from the split majority spin Fermi con¬ 
tour. We do not fully understand the origin of the weak 
peaks marked by open symbols in Fig. 3(b). They might 
stem from magnetic breakdown between the split con¬ 
tours [211 [22]; similar phenomenon has been invoked to 
explain anomalous SdH frequencies seen in bilayer elec¬ 
tron systems confined to double-QW samples in i?|| [H]. 

We summarize, in Fig. 3(c), the results of the Fermi 
contour calculations and the measured SdH frequencies 
(red and green squares), normalized to The cal¬ 

culated frequencies for majority- and minority-spin con¬ 
tours, which are equal to the calculated Fermi contour 
areas multiplied by /i/(27r^e), halve at i3||-values that 
mark the splitting of the respective contours. This jump 
in frequency reflects the fact that, for small H||, the plot¬ 
ted curves are based on the areas enclosed by the unbro¬ 
ken Fermi contours, whereas for larger H||, they are based 
on the areas of each of the split contours. There is good 
overall agreement between the measured and calculated 
SdH frequencies [23] . This suggests that SdH oscillations 
indeed show the spin-dependent splitting of Fermi con¬ 
tours, corroborating our interpretation of the two kinks 
observed in the H||-magnetoresistance (Fig. 2(a)). 

Having established the Fermi contour splitting through 
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FIG. 3. (color online) (a) SdH oscillations measured in the 
reference region of the Hall bar as B|| increases; traces are 
shifted vertically for clarity, (b) FT spectra of the SdH oscil¬ 
lations. The dotted and solid black lines show the expected 
positions of the spin-unresolved and spin-resolved FT peaks at 
Bjl =0, respectively. The FT signal to the left of the vertical 
lines indicated by ^10 and -h50 is affected by the Hamming 
window used in the Fourier analysis and is shown suppressed, 
(c) Summary of the FT peak positions normalized to fsdHj 
the frequency at i3|| = 0. Closed squares represent the mea¬ 
sured frequencies. The frequencies predicted by the calcula¬ 
tions for the spin subbands are shown as green and red lines. 
(A)-(E) mark the H|| values from Fig. 1 calculations. 



i3| I-magnetoresistance and SdH oscillations, we now turn 
to COs data measured in the modulated regions of the 
Hall bar. The magnetoresistance trace of Fig. 4(a), taken 
as a function of purely B±, is representative of such COs 
exhibiting pronounced minima at the electrostatic com- 
mensurability condition 2Rcfa = i — 1/4 [T^l [HHH], 
where i = 1, 2, 3..., Rc = kp/eB^ is the real-space 
cyclotron diameter, and a is the period of the poten¬ 
tial modulation [kp is the Fermi wave-vector perpendic¬ 
ular to the current direction). The frequency of COs, 
fco = 2hkp/ea, directly measures kp. Note that the 


FIG. 4. (color online) (a) GOs along [110] from the pat¬ 
terned region. Vertical lines mark the expected positions of 
the GOs resistance minima. The corresponding trace from 
the unpatterned reference region clearly shows no COs. (b), 
(c) COs from the patterned regions of the L-shaped Hall bar 
along [110] and [110] for different values of H||. In order to do 
Fourier analysis, we extract the oscillatory part of these traces 
as a function of 1/Bp by subtracting the background resis¬ 
tance. Panel (b) also includes three examples of this process. 
Traces are shifted vertically for clarity, (d), (e) Normalized 
FT spectra of the COs data shown in (a) and (b), respectively. 
The vertical dotted lines mark /qq (see text). 
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FIG. 5. (color online) Summary of the Fermi wave vectors 
(kp) deduced from the positions of the COs’ FT spectra for 
the two Hall bar arms. Blue and red symbols represent the 
experimental data for and respectively. (A)- 

(E) mark the B|| values from the calculations of Fig. 1. 

very high mobility of our sample leads to a large number 
of oscillations, up to f > 12. 

The magnetoresistance data from the [110] and [110] 
Hall bar arms are shown in Figs. 4(b) and (c). In each 
figure, the bottom traces, taken in the absence of i?||, 
exhibit high quality COs. As i?|| is increased, there is 
an obvious change in the periodicity of COs which is 
better seen in the FT spectra of Figs. 4(d) and (e). The 
bottom FT spectrum from each of these figures exhibits a 
single peak whose position (~ 0.35 T) is consistent with 
the commensurability frequency f^Q = 2hkF/ea = 0.35 
T (dotted line) expected for a circular, spin-degenerate, 
Fermi contour with kp = yj2'Kn. With increasing i?||, 
this peak / moves to higher frequencies in the FTs for the 
[110] Hall bar arm (Fig. 4(d)) and to lower frequencies 
in the [110] arm (Fig. 4(e)), suggesting that the Fermi 
contour is getting elongated. However, at = 12.5 T, 
a new peak / 1/2 emerges at approximately //2 in Fig. 
4(d). This indicates that the elongated contour has split 
into two smaller ones. As is increased farther, / 1/2 
develops to be the strongest feature of the FT spectra 
in Fig. 4(d). In contrast, / becomes progressively less 
pronounced and vanishes at i?|| = 16 T. 

Figure 5 summarizes the values oi kp extracted from 
the FT frequencies and also from the calculated Fermi 
contours for which we take the extrema along [110] and 
[110] (see the left inset) and plot them with bold red and 
blue lines, respectively. However, after the splitting, kp 
along [110], defined as shown in Fig. 5 right inset (thin 
blue line with an arrow), is represented by a thin blue 
line. We also plot half the length of the major-axis along 
[110] of the split contour (see the right inset) by a bold 
blue line. Qualitatively, the measured values oi kp show 


good agreement with the calculations, suggesting that 
the peak / 1/2 comes from the split Fermi contour. Cal¬ 
culations (see Fig. 1) also show that the extreme sizes of 
the contours for the two spin species always remain very 
similar. This explains why, unlike the SdH oscillations 
data, COs do not resolve the two spin Fermi contours 
m Another key point of Fig. 5 is that the elonga¬ 
tion of the Fermi contour deduced from the COs data 
is smaller than what the calculations predict. A similar 
discrepancy was previously observed in other 2D electron 
and hole systems [ini HU |30] . 

Another noteworthy feature of the COs data is that, 
even after the Fermi contour splits at H|| = 12 T, we 
appear to still follow / up to i?|| ~ 16 T (marked by 
solid blue circles in Figs. 4(d) and 5). To explain this, 
we propose a magnetic breakdown like scenario [mug 
where, even though the Fermi contour is split into two 
pieces, we still observe COs because a small portion of 
electrons jump between the split contours and complete 
the elongated orbit. In this context, magnetic breakdown 
in fc-space implies that, in real space, there is tunneling 
between layers which are separated because of strong B\\ 
(see Fig. 1(E) inset). 

We acknowledge support through the NSF (Grants 
DMR-1305691, ECCS-1508925 and MRSEC DMR- 
1420541), the Gordon and Betty Moore Foundation 
(Grant GBMF4420), and the Keck Foundation for sam¬ 
ple fabrication and characterization, the DOE BES 
(Grant DE-FG02-00-ER45841) for measurements and 
the NSE (Grant DMR-1310199) for calculations. This 
work was performed at the National High Magnetic Eield 
Laboratory, which is supported by NSF Cooperative 
Agreement DMR-1157490, by the State of Florida, and 
by the DOE. Work at Argonne was supported by DOE 
BES under Contract DE-AC02-06CH11357. We thank S. 
Hannahs, T. Murphy, A. Suslov, J. Park and G. Jones at 
NHMFL for valuable technical support. 


[1] N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Holt, Rinehart and Winston, Philadelphia, 1976), Chap¬ 
ter 12. 

[2] See, e.g., R. M. Potok, J. A. Folk, C. M. Marcus, and V. 
Umansky, Phys. Rev. Lett. 89, 266602 (2002). In trans¬ 
verse magnetic focusing experiments reported in this 
work, an anomalous shift of the focusing peaks’ positions 
is seen in Fig. 1 when H|| = 7 T. This shift, which is not 
discussed by Potok et al, stems from the i?||- induced 
distortion of the Fermi contour. 

[3] L. P. Rokhinson, V. Larkina, Y. B. Lyanda-Geller, L. N. 
Pfeiffer, and K. W. West, Phys. Rev. Lett. 93, 146601 
(2004). 

[4] J. P. Eisenstein, T. J. Gramila, L. N. Pfeiffer, and K. W. 
West, Phys. Rev. B 44, 6511 (1991). 

[5] A. Kurobe, 1. M. Castleton, E. H. Linfield, M. P. 
Grimshaw, K. M. Brown, D. A. Ritchie, M. Pepper, and 


















5 


G. A. C. Jones, Phys. Rev. B 50, 4889 (1994). 

[6] J. A. Simmons, S. K. Lyo, N. E. Harff, and J. F. Klem, 
Phys. Rev. Lett. 73 , 2256 (1994). 

[7] J. A. Simmons, N. E. Harff, and J. F. Klem, Phys. Rev. 
B 51, 11156 (1995). 

[8] N. E. Harff, J. A. Simmons, S. K. Lyo, J. F. Klem, G. S. 
Boebinger, L. N. Pfeiffer, and K. W. West, Phys. Rev. B 
Phys. Rev. B 55 , R13405(R) (1997). 

[9] M. A. Blount, J. A. Simmons, and S. K. Lyo, Phys. Rev. 
B 57 , 14882 (1998). 

[10] T. S. Lay, X. Ying, and M. Shayegan, Phys. Rev. B 52 , 
R5511(R) (1995). 

[11] T. Jungwirth, T. S. Lay, L. Smrcka, and M. Shayegan, 
Phys. Rev. B 56, 1029 (1997). 

[12] D. Weiss, K. von Klitzing, K. Ploog, and G. Weimann, 
Europhys. Lett. 8, 179 (1989). 

[13] E. Skuras, A. R. Long, 1. A. Larkin, J. H. Davies, and 
M. C. Holland, Appl. Phys. Lett. 70 , 871 (1997). 

[14] A. Endo, S. Katsumoto, and Y. lye, Phys. Rev. B 62 , 
16761 (2000). 

[15] D. Kamburov, M. Shayegan, R. Winkler, L. N. Pfeif¬ 
fer, K. W. West, and K. W. Baldwin, Phys. Rev. B 86, 
241302(R) (2012). 

[16] D. Kamburov, H. Shapourian, M. Shayegan, L. N. Pfeif¬ 
fer, K. W. West, K. W. Baldwin, and R. Winkler, Phys. 
Rev. B 85 , 121305(R) (2012). 

[17] D. Kamburov, M. A. Mueed, M. Shayegan, L. N. Pfeiffer, 
K. W. West, K. W. Baldwin, J. J. D. Lee, and R. Winkler, 
Phys. Rev. B 88, 125435 (2013). 

[18] R. Winkler, Spin-Orbit Coupling Effects in Two Di¬ 
mensional Electron and Hole Systems (Springer, Berlin, 


2003). 

[19] The calculations for = 0 predict a slight (<~ 5%) oc¬ 
cupation of the upper subband, which we do not observe 
experimentally. 

[20] The 7.08 T peak is observed because the By-range used 
in the Fourier analysis (0.2 < B±_ < 0.9 T) includes spin- 
resolved SdH oscillations starting at B± > 0.5 T. 

[21] N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Holt, Rinehart and Winston, Philadelphia, 1976), p. 
774. 

[22] J. Hu and A. H. MacDonald, Phys. Rev. B 46, 12554 
(1992). 

[23] Similar to a previous study |17| . the calculations, which 
neglect the exchange-correlation term between the spin 
species, predict a smaller spin-splitting than observed ex¬ 
perimentally. 

[24] R. W. Winkler, J.P. Kotthaus, and K. Ploog, Phys. Rev. 
Lett. 62, 1177 (1989). 

[25] R. R. Gerhardts, D. Weiss, and K. von Klitzing, Phys. 
Rev. Lett. 62, 1173 (1989). 

[26] C. W. J. Beenakker, Phys. Rev. Lett. 62, 2020 (1989). 

[27] P. H. Beton, E. S. Alves, P. G. Main, L. Eaves, M. W. 
Dellow, M. Henini, O. H. Hughes, S. P. Beaumont, and 
C. D. W. Wilkinson, Phys. Rev. B 42, 9229 (1990). 

[28] F. M. Peeters and P. Vasilopoulos, Phys. Rev. B 46, 4667 
(1992). 

[29] A. D. Mirlin and P. Wolfie, Phys. Rev. B 58 , 12986 
(1998). 

[30] L. Smrcka, arXiv: 1502.07129 


